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ABSTRACT 


We investigate the propagation of light rays in a clumpy universe constructed 
by cosmological version of the post-Newtonian approximation. We show that the 
linear approximation to the propagation equations is valid in the region z & 1 
even if the density contrast is much larger than unity. Based on a general order- 
of-magnitude statistical consideration, we argue that the linear approximation is 
still valid for z 'it 1. Then we give a general formula for the distance- redshift 
relation in a clumpy universe and derive an explicit expression for a simplified 
situation in which the effect of the gravitational potential of inhomogeneities 
dominates. In the light of the derived relation we discuss the validity of the Dyer- 
Roeder distance. Furthermore, we consider a simple model of an inhomogeneous 
universe and investigate statistical properties of light rays. We find that the result 
of this specific example also supports the validity of the linear approximation. 
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1. Introduction 


In the last several years large-scale structures of the universe has become 
an active area of research in cosmology. In spite of extensive theoretical as 
well as observational efforts, we have not yet understood the formation of the 
structures. Further observational information such as more complete survey of 
galaxies at high redshifts is definitely necessary in order to improve the situation. 
Fortunately it is expected that rapid progress in the observational techniques and 
the appearance of new telescopes of the next generation will bring about vital 
information on the structure of the universe in near future and open a new era 
for the observational cosmology. 

On the other hand it seems that careful attention has not been paid for theo- 
retical aspects of the observational cosmology. One of the difficulties is that light 
rays from distant galaxies might have propagated through intergalactic space in 
which the density is much lower than the average density of the universe. Thus it 
is not clear at all if the averaged homogeneous, isotropic Friedmann-Robertson- 
Walker (FRW) metric may be used as an appropriate metric on which light 
propagates or not. In fact the averaged FRW metric coincides nowhere with the 
real inhomogeneous metric. Since light feels local metric not the averaged metric, 
there is no justification for using the FRW metric to calculate the propagation. 
Nonetheless it has been customary to compare observations of distant galaxies 
with the predictions of a FRW universe. Otherwise rather crude descriptions 
of inhomogeneities such as the Dyer-Roeder model (1, or Swiss cheese model are 
used to interpret the observational data. It is not known in what sense such 
descriptions approximate inhomogeneities of the real universe. 

There is also another problem. Even if one has a realistic description of inho- 
mogeneities, it is not straightforward to relate theoretical quantities calculated 
on the inhomogeneous metric with actual observables. Without such relations, 
we will not have the correct interpretation of observational data. In view of 
rapid progress in the observational side it seems urgent to develop a consistent 
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theoretical framework of observational cosmology. 

Recently a consistent theoretical derivation of important relations in the ob- 
servational cosmology such as the magnitude-redshift relation was given in a 
linearly perturbed FRW uni verse. IJ ' J) Also an approximation method for an inho- 
mogeneous universe beyond linearized theory which is applicable to the present 
clumpy universe as well as to a linearly perturbed FRW universe is developed! 41 
The purpose of the present paper is to study the light propagation in a re- 
alistic inhomogeneous universe constructed by the approximation method men- 
tioned above and to extend the derivation of the magnitude-redshift (or distance- 
redshift) relation in the linearized case to the nonlinear case in which the density 
contrast is much larger than unity. 

The paper is organized as follows. In §2, we review the method for con- 
structing an approximate metric for an inhomogeneous universe and clarify the 
condition for the cosmological version of the Newtonian approximation. In §3, 
assuming that the condition for the Newtonian approximation is satisfied, we 
consider the propagation of light rays in an inhomogeneous universe and argue 
that the linear perturbation can be applied to the light propagation equations 
even in a highly inhomogeneous universe. In §4, we give the basic formula for 
the distance-redshift relation in an inhomogeneous universe, derive an explicit 
expression for a simplified case and compare it with the Dyer-Roeder distance. 
Our distance coincides with Dyer-Roeder’s in the region z 1. Then we con- 
sider a simple model of an inhomogeneous universe and compute the probability 
distribution of fluctuations in the distance-redshift relation. We find the result 
is in agreement with the general discussion of §3 and supports the validity of 
the linear approximation, even for very high redshifts. Finally, §5 is devoted to 
conclusions. 
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2. Approximation of inhomogeneous metric 

2.1. General scheme 

In order to make this paper self-contained, we shall briefly explain the ap- 
proximation method for constructing the metric of an inhomogeneous universe in 
general relativity developed by one of us. W Since the calculational detail has been 
presented in Ref. [4], we discuss mainly the physical idea behind the method and 
present the results. The method applies for the nonlinear stage as well as linear 
stage as far as the metric deviation from the FRW background metric. As shown 
below, this does not of course impose the smallness of the density contrast. 

We assume that the spacetime considered here may be parametrized by two 
independent small parameters e and k. The e is associated with the amplitude of 
the gravitational potential ((f)) generated by inhomogeneous distribution of mat- 
ter, <f> ~ e 2 . The k is the ratio between the typical scale of the inhomogeneities 
(f) and the scale of the background spacetime (L), it = t/L. The relative size 
of e and k depends on the system we have in mind. Since the metric fluctuation 
is generated by the density fluctuation Sp via Poisson equation, the density con- 
trast 6p/pb may be evaluated from A<j>/(Gpi) ~ e 2 //c 2 , where pi is the averaged 
density (see below). Thus the linear and nonlinear stage may be characterized 
by the condition k » e and e^> k, respectively. For example, if we take a super- 
cluster whose size is about 30 Mpch -1 , then k will be ~ 30/3000 ~ 10 -2 , where 
L ~ 3000 Mpch -1 is the present horizon size. We do not know the order of the 
gravitational potential for such a system, but the density contrast seems to be 
of the order of unity. Thus the gravitational potential would be e 2 ~ k 2 ~ 10~ 4 
for such a system. The size of supercluster seems to be the boundary between 
the linear and nonlinear regions. We also note that typical values of e and «c for 
galaxies are e ~ 10 -3 and k ~ 10 -4,5 . 

We make the following ansatz for the metric, 

Qtiv = > ( 2 * 1 ) 
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where a is the scale factor which describes the averaged global expansion and 
is assumed to be a function of the conformal time rj. It is also assumed that 
o'/a = 0(1 /L) where the prime means the derivative with respect to rj. The 
h’s are supposed to be generated by inhomogeneous distribution of matter and 
by possibly gravitational waves. We do not consider the latter possibility in this 
paper and assume that hp U = 0(e 2 ) and hp UiP = 0(e 2 /£). We assume that the 
spacetime considered here reduces to the closed, flat or open FRW spacetime 
depending on the curvature of the spatial section K — +1,0 or — 1, respectively, 
when the matter distributes homogeneously and hp U vanish identically. Thus the 
7p„ is the standard metric for one of the FRW universes. 

The above ansatz for the metric is used to expand the Einstein equations in 
terms of e and k as follows: 

3{2(^) + K}(K’ n + h) + 3{(^)’ + *> 

a 2 a ( 2 . 2 ) 

f 2 j j • 

{3(— ) + 2K}h* + (-)( S” + h) ' + , (2.3) 

a a Z Z ip 

(2— - (-) }(h^ + 7 ij h™ - 7 ij h) + K(& - W’h) 

a a 2 

+ {(-) +2--K}^ + -(2h^-h ii ^+ly ii h^) ( 2 - 4 ) 

a a a 2 

where we have used the trace reversed metric perturbation defined by h = 
h fU ' — j7 and worked in the harmonic gauge = 0. The indices on h 

and h are shifted by the background metric 7^ and the bar indicates a covariant 
derivative with respect to jp U . = a i T fiV +t ,iV is the total effective stress energy 
pseudotensor. The i* 1 " consists of terms quadratic in h and may be interpreted 
as a gravitational stress energy pseudotensor. In deriving the above equations, 
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we have neglected terms of order higher than 0 (c 4 /k 2 ) with the assumption 
that e 2 4Z. k which we assume throughout the present paper. Provided that 
one focuses on cosmological problems in which structures of interest are above 
galactic scales, this assumption gives practically no restriction for applicability 
of the present scheme to the real universe. 

The equations for the averaged global expansion is obtained by taking the 
spatial average of the above equations. The results may be written as follows: 

A’ = nr < 2 - 6 ) 

7 “TT (2 ' 6) 

+ (-4 - 2 (i) ! - 2jr ) <**> = i6 ” a ■ < 2 - 7 ) 

In deriving the above equations we have assumed ( r ’**) =0 expressing no coherent 
motion over the volume to be averaged and we have required that {h 1 ™) = (h%) = 
0 by choosing appropriate time variable and the scale factor. The spatial average 
of the line element takes then the following form, 

(ds 2 ) = a 2 {-dr? + ( 7 ii + (hiMdx'dxj] , (2.8) 

Thus (hij) express the deviation from the isotropic expansion due to the inho- 
mogeneities (t* j ) and the averaged spacetime expands anisotropically except if 
(hij) vanishes identically. Equations (2.5) and (2.6) axe the same with the equa- 
tions of the FRW model except that the source terms are replaced by the total 
effective stress energy pseudotensor including gravitational contribution. Thus 
the effect of local inhomogeneity on the global expansion may be expressed by 
the effective density p e ff — a 2 (r’ w ) and the effective pressure p e ff = jo 2 (r£). 
The equations which determine the local metric may be derived by subtracting 
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the above averaged equation from the original equations (2.2) ~ (2.4). These are 
as follows: 

Uh™ = -16irG(r im - (r™)) 

+ 7(*j? - *i|,) + 3 ( 2 ( 7 / - K)(h-» + S‘) , 

□ £* = -16t tGt* + ^•(h’wl* + h^) + 6(^) h * , 

~ ij a! - 

□ fc = — 167rGf ,J H 2^i„ 

a 1 ' 

2 

+ 4(^W + iy^J) - 2{(— ) - 2- - Jf>r , 

3 a a 

□ s* = -16<rG(T t * - (rf)) 

+ 7(^1, - <) + <{ (7) 2 - + *(«" - At) , 

where % = — j 7*^* I s the spacial trace free part of the perturbation and 

may be regarded as the gravitational wave degrees of freedom. Above equations 
are of course supplemented by the equations of motion = 0 or 

+ ^C 62 "" 1 + T W T) + (^ pk + i/ipj' 1 ) 2 ’'" “ (M~ + ^ lM r) = 0 . (2.13) 

Then the calculated perturbations hf iv are used to construct the local metric 

da 2 = a 2 (? 7)(7 mv + VO , (2.14) 

where a and ( h p V ) are calculated from the global equations (2.5) (2.7). 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


i 
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2.2. Linear and post- Newtonian approximation 


In this subsection, we construct the linear and post-Newtonian approximation 
from the equations derived in the above. First we give a characterization for the 
linear and nonlinear regimes. Since K 1 ™ ~ e 2 and Gpb ~ l/£ 2 , pj the averaged 
density, by definition, we can easily obtain the order of the density contrast from 
the local equation (2.9), 


6p e 2 
Pb 


(2.15) 


Thus the linear and nonlinear regimes are characterized by the condition k e 
and e ^ k (>> e 2 by assumption), respectively. 


We first consider the linear approximation which applies for the linear regime 
k ^ e. In this regime the dynamical timescale of the density fluctuations will be 
of the order of the timescale of cosmic expansion. Then the equations of motion 
is used to evaluate the order of the velocity as e 2 /«. Thus we have the following 
ordering in the linear regime, 


Sp e 2 _ , e 2 


If" ~ e 2 < 1, — ~ -x < 1, v* < 1 . 


Pb 


K* 


K 


(2.16) 


Every perturbed quantities are much less than unity and their second order terms 
are safely ignored. For example, t**" may be approximated by a i T fil '. This is of 
corse what we call the linear approximation. If we take the perfect fluid form 
for the stress energy tensor, we may neglect the spacial trace free part because 
these are second order in velocity and thus the averaged spacial trace free part of 
the metric perturbation ^ may be neglected. Thus the spacetime expands 
isotropically in this case. The expansion equations (2.5) and (2.6) reduce to the 
usual equations for a FRW model. 

Next we consider a cosmological version of the post-Newtonian approximation 
which applies for the nonlinear regime e k. In this case the dynamics of the 
density fluctuation is totally governed by its self gravity and thus the dynamical 
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time scale will be the Newtonian time defined by r = eql* 1 The equations of 
motion then give the proper relation <f> ~ v 2 where <f> is the Newtonian potential 
generated by the density fluctuation. Thus we have the following ordering in the 
nonlinear regime, 

<t> ~ e 2 < 1, h ij ~ e 4 < 1, v* ~ e < 1 , — ~ ^ > 1 . (2.17) 

p b K 2 

This ordering allows us to neglect terms like (a 1 /a) 2 h ttv in the equations (2.9) 
~ (2.11) and to solve the equations perturbatively. The lowest order equations 
are just the Newtonian equations and the next order is order e 2 smaller than the 
Newtonian order. Thus we may safely use the Newtonian approximation as far as 
e is sufficiently small, but much larger than ic. In such a region the approximated 
line element is given by 

ds 2 = a 2 {rji) [-(1 + 2<b)dri 2 + (1 - 2^)7 ijdx'dx*] , (2.18) 

where a and ¥ axe determined by the following equations in the lowest order: 

(— )* = - K (2.19) 

CL O 


A * = 4tt Ga 2 6p , (2.20) 

(») 

where A is the Laplacian operator in the spacial section. Thus is the Newtonian 
potential generated by the density inhomogeneity. The equations of motion in 
the lowest order are the usual conservation law and the Newtonian equation: 

p' +3-p + (pv i ), i = 0, (2.21) 

a 



( 2 . 22 ) 


where v l = dx % /dr\. In particular, Eq.(2.21) implies Pb oc a~ 3 at this order. Thus 
the conventional FRW model is justified as a model for the averaged background 
universe. The higher order corrections are calculated elsewhere. 
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It should be pointed out that we may write the line element in the above 
form (2.18) even in the linear approximation by choosing an appropriate gauge, 
namely the so-called Newtonian gauge! 2 '* 1 It is known in the linearized theory 
that the potential satisfies the same Poisson equation (2.20) and is interpreted 
as the Newtonian potential. The difference between the linear approximation 
and Newtonian approximation is then the equations of motion. We may thus use 
the above expression for the line element in the linear as well as the nonlinear 
regions. 


3. Light propagation in an inhomogeneous universe 

We now consider the propagation of light rays in an inhomogeneous universe 
constructed by the previous method. Thus the line element takes the form (2.18) 
in the linear as well as the nonlinear regimes. We may safely adopt the geometric 
optics approximation in the cosmological context. Since the constructed metric 
takes the form g M „ = a 2 ( 7 M „+h #4 „) = a 2 £ M „ and the light propagation is unaffected 
by conformal transformations, it is rather convenient to work in the conformally 
related spacetime g^ v . In the following, quantities with tilde axe quantities in the 
conformally related world. 

The basic equations for the propagation of light rays are 

^" + ^“^=0, (3.1) 

= -R a 0k a hP - - 2*’ , (3.2) 

= ~C a p^ a h p k <T — O&af) ) (®*®) 

where the first equation is the null geodesic equation and the second and third are 
the equations for the expansion 6 and shear tensor a a p of the bundle of light rays. 
Our strategy here is to show first that fluctuations in Eqs.(3.1) and (3.3) due to 
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the presence of inhomogeneities axe small. Then we consider Eq.(3.2) and look 
for conditions for the validity of the linear approximation for the propagation of 
light rays. 

We write the 4 -momentum of the light rays as follows: 

k p = k* + 6k * (3.4) 

where k p is the 4-momentum in the unperturbed metric = 0 and satisfies the 
geodesic equation, 

+ T^k'k* = 0 (3.5) 

where Tpo is the Christoffel symbol of the metric 7 ^. This equation then tells 
us that k ** = 0(1). The 6k p satisfies the following equation, 

^r6k» + + 2T p <r k f> 6k 9 = 0 , (3.6) 

CL A 

where we have neglected higher order terms and 6Tp<x is given by 

+ < 3 - 7 ) 


From this expression we find that — 0(e 2 /£). 

As for the shear, it vanishes on the unperturbed background since the Weyl 
tensor is identically zero for 7 ^, provided the light rays are originally shear-free. 
Hence the propagation equation at the lowest non-vanishing order is 

= ~C a pp tr k p k <r — 9<r a p , ( 3 - 8 ) 

where we have C ap p a = 0 (e 2 /£ 2 ). 
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Now consider the propagation equation for the expansion. The Ricci tensor 
is approximately written as 


Rfitz — Rfii* + & Rfi 1 / , ( 3 . 9 ) 

where R^ = 2K'yij6 t fl 8i is the Ricci tensor of the unperturbed spacetime and 
SR^ is given by 


SRfHf ^pa\iiu) • (3.10) 

Note that 6R ltt/ = 0{t 2 /l 2 ). We decompose the expansion as 

9 — 9 -t- 69 , (3.11) 

where 6 is the expansion on the unperturbed background. Then 6 and 86 satisfy 
the following equations, 

= - 2 * - . ( 312 ) 

^ 69 = —6{R tlv k ,l k ,, )\ - 969 - 1 69 2 - 2a 2 , (3.13) 

where 

6(R fiV k' i k u ) x = 6R liV k> t k v + *k4?{ kitx') , (3.14) 

with 6x * being the perturbation of the geodesic path and we have used the fact 
that the shear vanishes in the unperturbed FRW universe and assumed that the 
perturbation 6 k* is small, which we shall verify shortly. Note that the expansion 
on the unperturbed background is 9 = 0(1/A). 

First we consider 8k 11 . Imagine a universe filled with objects (galaxies) of 
size l whose density is 8p = (e 2 //c 2 )pj. Hence the mean separation distance is 
r 0 = (e/«) 2 / 3 ^. Then for a geodesic affine distance of A, the light gravitationally 
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encounters such objects N g = A/ro = e~ 2 ^K~ 1 ^(X/L) times in average. In each 
encounter, the integral of the geodesic equation gives the contribution (e 2 /6) x 
b = e 2 to where b is the impact parameter (6 ^ ro). However, since the 
sign of each contribution will be random, the total contribution to 6k ** will be 
y/Ngt 2 = e 4 / 3 (e 2 //c) 1 / 6 (A/£) 1 / 2 , which is always much smaller than unity, since 
t 2 < s by assumption and A L (note that the affine distance to the source 
object A is bounded from above due to the fact that the age of our universe is 
finite). It should be mentioned that since the contribution of each encounter is 
independent of b, the direct encounters which are much rarer events are totally 
unimportant in this case. Thus we conclude that the linear approximation is 
valid for the evaluation of 6k ft . 

As for the magnitude of the shear, the same argument as above leads to the 
estimate a 2 ~ (c 2 /x)(A/i 3 ) as the contribution from gravitational scattering, 
where we have assumed that the mean impact parameter (b) is of the order ro. 
The contribution from the direct encounters can be similarly estimated by noting 
that the average number of encounters is Nj = (£ 2 /rj *)A = («/e 2 )(A/£) with each 
encounter contributing (e 2 /£ 2 ) x £ = e 2 /i with random sign. The result turns out 
to be the same as that of gravitational distant encounters. Hence we conclude 
that a 2 ~ (e 2 //c)(A/Zr 3 ). 

Finally let us consider the expansion. From the above estimate for the shear 
and from the fact 0 = 0(1/A), we find 66/6 = (e 2 //c)(A/£) 3 < 1 as the con- 
tribution of the shear. Thus the linear approximation for the evaluations of the 
shear and its contribution to 60 is justified. To evaluate the contribution from 
the perturbed Ricci tensor, it is important to know the explicit form of 6R lu/ k ,i k u 
(see Eq.(4.3) below); 

«p) • (3-15) 

Hence from the Poisson equation (2.20), the dominant contribution from 6R tlv 
to 60 comes only from regions where 6pj p is non-vanishing, i.e., there is no 


6R ltv k ll k v = 2 A $ + 


°( 
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contribution from distant gravitational encounters at the leading order. 

The contribution from the direct encounters with galaxies can be estimated 
similarly as above. We obtain 60 ~ (t 2 /l 2 ) x l x JV^, where note that each 
contribution to 66 is negative definite (since 6p/pi is positive definite) in this 
case, hence adds up secularly as oc Nj. This gives 60 = 0{X/L 2 ) and we obtain 
60/0 = 0(\ 2 /L 2 ). Hence the linear approximation is valid if (A /L) 2 1, i.e., 

z 2 <£. 1. On the other hand, we have = 0(1/L 2 ) in intergalactic space. 

Hence it gives 60 ~ (1/L 2 ) x A which is the same as the contribution from galaxies 
in magnitude but has the opposite sign (since 6p/ p\, is negative definite). There- 
fore, provided the number of encounters of a light ray with galaxies is sufficiently 
large (N& ~ 30(A /L) for galactic scale objects), we can expect these two contri- 
butions to cancel each other on average. The resudial fluctuation in 60 will be 
proportional to y/Nl > M ~ (e 2 /l 2 ) xlx y/Nj, hence (60 /6) 2 ~ (e 2 / /c)(A/ L) 3 1. 

Note that the assumption e 2 « guarantees Nj 1 for z ^ 1. As a result the 
linear approximation will be valid even for A = O(L), i.e., for :>1. 

The above argument may sound too naive. However, as for the shear con- 
tribution, there exists a more detailed theoretical argument [T1 which supports our 
result. As for the perturbed Ricci tensor contribution, in the next section, we 
shall investigate statistical properties of light rays in a simple but reasonable 
model universe and show that the results are indeed consistent with the above 
argument. Further, recent numerical calculations by Watanabe 1 * 1 also seem to 
support the validity of the linear approximation. To summerize, we conclude 
that except for a statistically very rare kind of light rays, the linear approxima- 
tion can be safely used to study the propagation of light rays in an inhomogeneous 
universe in which the density contrast is much larger than unity as long as one 
focuses on a region sufficiently smaller than the horizon scale, i.e., z < 1, and it 
is a very good (if not the best) approximation even for z >■ 1, provided that our 
approximation based on the expansion in terms of e is applicable to the universe, 
i.e., e 2 k. 
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4. Distance-redshift relation 


4.1. Basic formulas 


As discussed in §2, the approximate metric (2.18) gives a sufficiently accurate 
description of the real universe. Further, we have seen in §3 that the linear 
approximation can be used to investigate the propagation of light rays in the 
region z & 1 (and probably even for z £ 1, provided one is interested in average 
light rays). These imply that most of the arguments given up to §4 of Ref. [2] 
hold also for a highly inhomogeneous universe, under the assumption that e 2 k, 
since only geometrical considerations but non of the Einstein equations were used 
there. The only modifications we have to make are to include the contribution of 
the shear in the propagation equation for the expansion and to retain terms of 
0(v 2 ) in the expressions for perturbed four- velocities of the source and observer. 
Note that 0($) ~ 0(v 2 ) in the Newtonian situation. 

It follows that the distance-redshift relation in an inhomogeneous universe, 
written in the Newtonian gauge, is expressed as 12 ’* 1 


Sdijz, 7 ») _ ^( 3 , 7 *) 
d A {z) 


= coth V—K^t + I , 


(4.1) 


with S\ a and I given by 


/ / „t \ -1 \ 


1 

l ^ 

f 

a, = ((a). ~ A T 

l 

¥ - v»7* - -« 3 
2 

+ 2 /<iAtf'(A) 

0 i 
0 


+ 2 



¥ - $0 + 



+ A a (¥ + v*7* + -v 2 )o > 


(4.2) 
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I =— != J dX sinh 2 V^KX (coth V^KX - coth V^KX,) 
0 

xji «(*„*<■*•'),+*»} ; 


\s(R^k“k-) x =A 4 - (4" + 2^4') 


+ 2ir( 


A. 

(4_„ j y_I^), + 2 JdXt 4*(A,) 


( 4 . 3 ) 


where Sdi/di (Sd^/dji) is the anisotropy in the luminosity (angular diameter) 
distance to a given redshift z measured in the direction of 7’, v 2 = r,v‘, X t = 
r](z) — r/o is the conformal distance to the redshift z in the averaged background 
universe, the suffix s denotes a quantity at the source of redshift z, the suffix 0 
for z — 0, ^'(A) = d'S (77(A), x t (X))/dri, etc.. Note that di = (1 + z) 2 ^ holds 
for arbitrary spacetimes J 11 01 ■ which guarantees tfdfc/di = Sd/^fd^ (= Sd/d). The 
difference between the case of linear density perturbations and the non-linear case 
is that and v * are determined from the Newtonian equations (2.21) and (2.22) 

in the latter, while they are determined from the linearized Einstein equations in 
the former. 


For distances in the range ^ < A t & L, the examination of the order of 
magnitude of each term in Eq.( 4 . 1 ) ~ ( 4 . 3 ), in the non-linear cue, shows that 
the leading terms, up to the order of e 2 /#e, are those involving A ’® r , £ 2 &nd Vij*. 
Neglecting the other terms we find 


Sd 

d 


- V'-tf ®oth y/^KX, { (^7*), - (^7%} + (*•?*), 

A. 

— —== J dX sinh 2 y/—KX (coth y/—KX — coth y/—KX t ) 

0 

x {AirGSpa? + d 2 } , 


( 4 . 4 ) 
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where the Poisson equation (2.20) was used to replace A '® f by Sp. Further com- 
parison of the remaining terms, with the same strategy we took in the previous 
section, shows the first term (oc {(v»7*)j — (v«7*)o}) is generally most important 
for 


z ~ \,/L £ e 1 / 3 = z cr , 


(4.5) 


and the density term (oc Sp), otherwise. However, the other terms can be impor- 
tant at high redshifts if the integral of the density term gives zero in average as 
we discussed in §3. We note that the first term describes the usual Doppler effect 
due to the peculiar velocities of the source object and the observer, which is used 
in the determination of the peculiar velocity field on large scales (^ 50Mpch -1 ). 
Note also that with typical values of e for galaxies, one has z cr 10 *, or a 
distance of ~ 300Mpch -1 . 

Now let us discuss the relation of Eq.(4.4) to the Dyer-Roeder distance! 11 
To do so, we must assume the same physical situation as the one assumed in 
deriving the Dyer-Roeder distance. That is, we assume that light propagates 
only through the intergalactic space where the density is uniform and given by 
PIG = a Pb (0 < a < 1). Hence Sp = —(1 — a)pj. In addition we assume that 
the potential gradients can be neglected in the intergalactic space so that there 
is no shear contribution. Further we consider a spatially flat universe (K = 0) 
for simplicity. Then Eq.(4.4) is easily evaluated to yield 


Sd 

d 


2{y/l ~+~Z — 1 ) 
+ 3(1 - a) 


[(«*)» - (v»)o] 7' + (vi)*7‘ 


VTTT + 1 


L v i+2 — i 


ln(l + z) — 4 


(4.6) 


In accordance with the general argument used to derive Eq.(4.5), the above equa- 
tion clearly shows that the Doppler term dominates for z ^ n 1 / 3 ~ e 1 / 3 but 
becomes unimportant for z v 1 / 3 . 
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Let us concentrate on the case when the peculiar velocities are negligible and 
compare our result with the Dyer-Roeder distance. The Dyer-Roeder (luminos- 
ity) distance for the spatially flat universe is given by ul 

d DR (z; Q ) = ^(l + *) S/4 [(1 + zf 1 ' ~ (1 + z)-* 1 *] , (4.7) 

where /3 = -J25 — 24a. It is then easily shown that our result (4.6) coincides 
with Eq.(4.7) in either of the limit 1 - a < 1 or 2 < 1. In fact, the coincidence 
is much more impressive than one would formally expect; the relative error is 
less than 0.5% at z = 1 and about 10% at z = 5 even for the extreme case of 
a = 0. Of course, this coincidence is not accidental. If we examine the essential 
assumption which lead to Eq.(4.6)(with the velocity terms neglected), we find it is 
the fact that the perturbation in the affine parameter distance SX t is negligible, 
*.e., one can use the geodesic equation on the background FRW spacetime to 
relate the redshift of the source object with the affine distance. This is just the 
assumption used to derive the Dyer-Roeder distance. Hence apart from the fact 
that we employed the linear approximation in evaluating the expansion of light 
rays, there is no essential difference between the two. 

The above discussion, together with the considerations given in §3, shows the 
validity and the limitation of the use of the Dyer-Roeder distance. First of all, it 
is valid only if e 2 <icso that the contribution of the shear can be neglected. Then 
it is valid for 1 £ z ^ z 6 * ~ e 1 / 3 if light rays which reach us came only through 
intergalactic space. Since the fact that Sk** is small is always true even for z ^ 1, 
the Dyer-Roeder distance is, in the mathematical sense, probably valid also for 
z 1. However, because the condition e 3 4Z. k implies that the expected number 
of encounters of a light ray with galactic objects are large; Nj ~ /c/e 2 1 for 

2^1 (see §3), it becomes increasingly rare for a light ray to pass only through 
intergalactic space. Consequently, the Dyer-Roeder distance becomes physically 
almost useless for 2 ^ 1 . 


19 



4.2. Statistical consideration 


We now focus on the contribution of the density term, which is seemingly 
not only dominant but also non-linear, to the propagation of light rays and 
investigate their statistical properties in terms of the distance-redshift relation, 
under the linear approximation, in a simple model of an inhomogeneous universe 
(see Ref.[ll] for a similar analysis). In particular, our main intention is to show 
that the results are consistent with the linear approximation. Hence our basic 
equation is 

j =(x)- 2 , = £/" (a - a ' )aS *- (4 ' 8) 

o 

We consider a spatially flat universe whose average density is in which stellar 
objects (galaxies) of equal mass are placed on a cubic lattice of constant comoving 
size, and the rest of space (intergalactic space) has a uniform density apf,. The 
galaxies are assumed to be spherical and uniform in density with a constant 
proper radius 12*. Then, we employ the so-called Wigner-Seitz approximation, 
frequently used in nuclear physics, and replace each cube of the lattice by a sphere 
of comoving radius tq. With this approximation, the Einstein equations in each 
Wigner-Seitz cell become 



SttG o 

=-3-p* a , 


A ¥ =47tG(1 — a)pi>a 2 


rga 3 


Rl 


d(R, — or) — 1 


r < ro , 


(4.9) 


where 6{x) is the step function and the center of the cell is chosen to be the 
origin. 


We divide the integral in Eq.(4.8) into intervals [A<, A* + AA] (AA ~ ro A,); 

Aj'fAA 

m 

(4.10) 
ii 

In order to avoid inessential complications, we further replace the sphere of radius 


/ = £/<; Ii = M l jzM j S*. 
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tq by a cylinder of radius tq and height AA such that the volumes of the two 
are equal and assume that a light ray comes into the cylinder parallel to its axis. 
Thus we take A A = (4/3)ro. Then in each interval, if the light ray passes through 
the galaxy, on average it will give the integral the contribution, 


Ai+AA 


/ 

Ai 


(») 


! ( ar o) 3 


R* 


dX A * =4ttG(1 - a)p b a 2 ^- xjx2 

ii; 3 a 

16 nti \ 2 ( ar o) 2 

=y? r G( 1 - a)p b a'r 0 — . 


(4.11) 


This occurs with the probability p = R 2 /(aro) 2 . On the other hand, regardless 
of whether the ray hits the galaxy or not, there is always a contribution to the 
integral given by 


Aj+AA 

J dA A $ = — 47rG(l — a)p b a 2 x ^ro 

Ai 
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= - y7fG(l “ «)pto 2 ro . 


Thus we may rewrite Ii as I{ = c,xj where 


ci =8(1 — a) 


7*q A,(A; — A«) 
A, ( 770 -Ai) 2 ’ 


x i = 


I — °* r ) ~ * 

LPi 


Pi = 


2 

IwF’ 


(4.12) 


and the back ground equations (p b oc a 3 , a oc rj 2 ) have been used. Then the 
probability distribution of Xi is given by 


■P»(*») = Pi& (*» + (1 - Pi)S(xi + 1) , 

and the probability distribution of J = ^ ciXi is expressed as 

P(I)dI = J fpXiP.txi) « (j - £ dI ' 


(4.13) 


(4.14) 


We note that Eq.(4.13) implies (xi) = 0, which in turn implies ( I ) = 0, the result 
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in accordance with the argument given in §3. 

To evaluate the probability distribution function P(I), we consider the char- 
acteristic function of P(I). We find 


OO 

F(q) = J P(I)e^'dI 


-n 


oo 

J dxi P(xi)e^ CiXi 


L-oo 


=e ~ iq c * exp ^ln [l + Pi (e^ Ci ^ - 1)] j 


Using the fact pi 1, the above can be approximated as 


F(q) = e ’*£« Ci exp 
Hence, at q 1, we obtain 




F(s) w ’*£‘“ exp ^ (ic tq - 


— exp 

— exp 


vi, 

A* 2pi 




where from Eqs.(4.12), the function f(z) is given by 


\ ( y/T+x-l\ 


(4.15) 


(4.16) 


(4.17) 




(4.18) 


where Ro = aot*o is the present proper radius of the comoving cell, *.e., half the 
mean separation distance of galaxies and Hq = 100 hkms -1 Mpc -1 is the Hubble 
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constant. On the other hand, for q — ► oo, the exponential factor in the exponent 
of F(q) oscillates rapidly and vanishes on average. Hence we find 


F(q) «e *®S» e * 

= e -ff(*) c +*9^o(*) ? 


(4.19) 


where g(z) and Iq(z) are given by 




0 + 2 + 1 


\50 kpc 

Io W =3(l-«)[^__i 


flo 

ln(l + z) - 4| , 


(4.20) 


where Jo is just the Dyer-Roeder part of the distance-redshift relation given in 
Eq.(4.6). It is also clear that g(z) represents the “optical depth” associated with 
encounters of light rays with galaxies! 71 

Now combining Eqs.(4.17) and (4.19), and taking into account the facts that 
(1) = jP{I)dI = 1 and (I) = JIP(I)dI = 0, we deduce that 


p(i) =+ J ni)'-**! 

(/-W) 2 ' 

2/(3) J 


.(1 _ e -<»(*)) 


■v/27 7f(z) 


exp 


(4.21) 


where I\ is given by 

A(3) = - 5i ^ rT /oW. (4.22) 

The interpretation of the above distribution function is very easy. The term 
proportional to the delta function represents light rays which never encounter a 
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galaxy, hence giving the Dyer-Roeder distance for I . Of course, the appearance of 
the delta function is due to the simplification we adopted for our model and not 
real. For example, any contribution from the shear or the density fluctuations in 
intergalactic space will develop a finite dispersion. The other term of a Gaussian 
form in Eq.(4.21) represents those light rays which encounter galaxies at least 
once. Thus I\(z) and f(z) are the mean value and the dispersion, respectively, 
of I for such rays. 

What we have to show is that the probability distribution function (4.21) 
is consistent with the assumption used to derive it, i.e., the linear approxima- 
tion. The consistency demands (J 2 ) < 1. The expectation value (I 2 ) is easily 
calculated to be 

</ 2 > = ^hr IoM ‘ +(1 ~ ■ ( 4 - 23 > 

It is easy to see that the first term proportional to I 2 is always very small at any 
redshift. Incidentally, this implies I\ is very small at any redshift (an order-of- 
magnitude evaluation shows |/i| m a* = 0(e 4 //c 2 ) at z ~ e 2 //c). The smallness of 
the second term will be guaranteed if f(z) 1 for any z. As given in Eq.(4.18), 
for characteristic values of R* and Ro, we find f(z) is always smaller than unity 
(note that Eq.(4.18) is indeed in the form f(z) = 0(c 2 /k), in agreement with 
the order estimate given in §3). Thus not only the general order-of-magnitude 
argument of §3 but also the specific model consideration of this section support 
strongly the validity of the linear approximation for study of the propagation 
of light rays in a highly inhomogeneous universe, and hence the validity of the 
formula (4.1). 
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5. Conclusions 


We have carefully investigated the propagation of light rays in a realistic 
inhomogeneous universe. We have found that the linear approximation to the 
propagation equations is valid even in a universe with very high density contrast. 
Then we have derived a general expression for the distance- redshift relation (4.1). 
It should be stressed that the relation has been derived totally within the frame- 
work of general relativity without any ad hoc assumption. 

Based on the drived relation, we have been able to clarify the validity and 
the limitation of the use of the Dyer-Roeder distance; it is numerically valid only 
for 1 £ z ^ e 1 / 3 if the light rays which reach us have traveled only through inter- 
galactic space. In more general situations, one has to use Eq.(4.1) to construct 
the distance-redshift relation. The equation can be evaluated explicitly once the 
density and velocity distributions are given throughout the region of interest. 

We have given some plausible arguments that thus constructed relation is 
valid for any redshift provided e 2 « is satisfied and the condition gives prac- 
tically no restriction for its applicability to the real universe as fax as regions 
of interest are above galactic scales. Although the validity of Eq.(4.1) should be 
checked more rigorously (for example, by a careful numerical analysis which takes 
full account of non-linearity in the propagation equations), if it is indeed justi- 
fied, it will play a fundamental role in the correct interpretation of cosmological 
observations. 
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